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■ Abstract Chandra Galactic Center Survey detected ~ 800 X-ray point- 

I ' like sources in the 2° x 0.8° sky region around the Galactic Center. In this 

| H I paper, we study the spatial and luminosity distributions of these sources 

, according to their spectral properties. Fourteen bright sources detected are 

^ ■ used to fit jointly an absorbed power-law model, from which the power-law 

' photon index is determined to be ^2.5. Assuming that all other sources 

H ; 

^ have the same power-law form, the relation between hardness ratio and HI 

column density Nh is used to estimate the Nh values for all sources. Monte 
Carlo simulations show that these sources are more likely concentrated in 
the Galactic center region, rather than distributed throughout the Galactic 
disk. We also find that the luminosities of the sources are positively cor- 
related with their HI column densities, i.e. a more luminous source has a 
higher HI column density. From this relation, we suggest that the X-ray 
luminosity comes from the interaction between an isolated old neutron star 
and interstellar medium (mainly dense molecular clouds). Using the stan- 
dard Bondi accretion theory and the statistical information of molecular 
clouds in the Galactic center, we confirm this positive correlation and cal- 
culate the luminosity range in this scenario, which is consistent with the 
observation (10'^^ ^ 10"^^ ergs s^^). 
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1 INTRODUCTION 



The central regions of galaxies are usually crowded by many celestial bodies with different 
physical properties. The center of our Galaxy is an ideal laboratory to study X-ray 
sources and the related high energy astrophysical processes. It has been the observational 
target of most X-ray satellites, such as ROSAT (Snowden et al. 1997; Sidoli, Belloni & 
Meregetti 2001), ASCA (Sakano et al. 2002) and BeppoSAX (SidoU et al. 1999). In 
these past observations, many X-ray point-like sources have been detected. The ROSAT 
survey covered a 3° x 4° sky region around the Galactic Center and detected 107 X- 
ray point-like sources in 0.1-2.4 keV. In ASCA observations, 52 point-like sources were 
detected. BeppoSAX has also observed 16 X-ray point-like sources in the central region 
of our Galaxy. Chandra X-ray Observatory {CXO) has much better spatial resolution 
(0.5 arcsec) than previous satellites (e.g. ROSAT 2Q arcsec). It can produce high quality 
images of the Galactic center. Chandra Galactic Center Survey (hereafter GCS) was 
carried out in July, 2001 (Wang, Gotthelf & Lang 2002). It consists of 30 separate ACIS- 
I observations (Obs. ID 2267-2296, energy band 0.2-10.0 keV) with a total exposure time 
of 94.2 hours (~ 340 ks). The sky area covered is 2° x 0.8°, about 280 pc x 110 pc 
(taking the distance of Galactic center to be 8.0 kpc). About 800 point-like sources were 
detected. 

A variety of analyses to explore the nature of these discrete sources have been done. 
Pfahl, Rappaport & Podsiadlowski (2002) discussed the possibility of wind- accreting neu- 
tron star nature of Chandra GCS sources, while Belczynski & Taam (2004) introduced 
another possibility for low luminosity (10^^ ~ 10^^ ergs s~^) source nature, the Roche 
lobe overflow accretion systems. Belczynski & Taam (2004) also argued that accreting 
neutron star systems are not likely to be the nature of the major part of the sources. 
Multiwavelength observations have been proposed on these sources to further explore the 
nature of these sources. ChaMPlane [Chandra Multiwavelength Plane Survey, Grindlay 
et al. 2005) focuses on the multiwavelength observations on the Chandra X-ray sources in 
the Galactic plane and bulge. They employed X-ray and optical surveys, and a foUowup 
spectroscopy and an IR identification program. Their recent results (Laycock et al. 2005) 
showed that high mass X-ray binaries (HXMBs) are not the dominant population (fewer 
than 10%) of the X-ray sources in the Galactic center. Bandyopadhyay et al. (2005) took 
an infrared survey on parts of the Chandra GCS region. Analysis on potential coun- 
terparts indicated that the sources may be accreting X-ray binaries. Recently, Liu & Li 
(2005) carried out an evolutionary population synthesis study on discrete sources in the 
Galactic center region. Their results showed neutron star low mass X-ray binaries could 
not account for most sources in Wang et al. (2002) survey. Muno et al. (2006) reanal- 
ysed the data of Chandra GCS and suggested that major part of these sources could be 
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Fig. 1 The total chi-squared val- 
ues for different photon indices. When 
r=2.5, the reduced chi-squared value 
is most close to 1. 
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Fig. 2 Fitted spectrum of CXO 
J174417.2-293945: Nh = (0.5476 ± 
0.0876) X 10^2 cm-2, r=2.5. 



CVs since the local Galactic neighborhood X-ray sources showed a CV population in the 
10^2 — 10'^^ ergs luminosity range. 

Despite of extensive work done previously, the nature of X-ray sources in the Galactic 
center is still not understood completely. Our work presented here is a statistical analysis 
of the sources in Wang's Chandra GCS, including their spatial distribution and flux 
characteristics. In §2, data analysis method and results arc described in detail, including 
spectral fitting and parameter determination. In §3, the spatial distribution of the sources 
is discussed. In §4, another potential nature of GCS X-ray sources is proposed. §5 is a 
summary. 



2 SPECTRAL FITTING & Nh DETERMINATION 

Chandra has excellent angular resolution. Although many GCS sources have only a few 
counts, the background counts arc even much lower. So for most sources, the detec- 
tion probability is close to unity. The efficiency correction is ignored here because the 
luminosity function of these sources is not discussed in this paper. 

Psextract procedure in the CIAO package is used to get the spectral files of the 
sources. The regions we used to extract energy spectra are circles chosen to contain 
all source photons, while minimizing the effect of nearby sources and background. The 
radii of these circles range from 6 pixels to 24 pixels. The background region is made 
of several circles of similar radii for these source. The spectra are re-binned to make 
each energy channel have no less than 25 counts for most sources. For other sources with 
smaller counts, each energy channel has no less than 10 counts. Then the Xspec in Heasoft 
package is used to fit the spectra with the absorbed power-law model, i.e. "phabs(po)" 
model in Xspec. There are three parameters in this model — the photon spectral index 
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r, the neutral hydrogen column density Nh and the normalization factor A. Because 
of the short effective exposure time for each field, the total counts for sources without 
significant pileup are not sufhcient for a high quality spectral fitting, in order to determine 
Nh and T simultaneously for each source; many other sources arc also too faint for any 
meaningful spectral fitting individually. We thus choose 14 sources, each with total counts 
of more than 80 and without significant pileup (the maxim pileup percentage of the 14 
sources is 8%). We assume that all sources have a similar value of F, thus attributing 
the different hardness ratios of these sources to different values of Nh- Our goal is to 
search for one value of F which best describes all 14 sources. Spectral fittings of these 14 
sources show their power-law photon indices ^2 — 3, indicating soft spectral properties. 
F values are tried from 1.5 to 3.5 with step of 0.2. For each trial value of F, i.e., we fit the 
absorbed power-law model by fixing F, the fitting generates one chi-squared value. Then 
the total chi-squared values of the 14 sources for all different values of F are calculated. 
We find that F = 2.5 corresponds to the minimum chi-squared value (Fig. 1), with a 
reduced chi-squared value of about 1.180 for 185 degrees of freedom. We thus conclude 
that F = 2.51q'5 f'^'' ^ 68.3% confidence interval. We believe this single photon index 
can describe most of the GCS sources adequately and no obvious selection effects exist, 
because these 14 sources distribute through the whole ranges of A^^f distribution acquired 
below (Fig. 4). Table 1 shows the spectral fitting results for these 14 sources with the 
photon index frozen at 2.5. 



Table 1 Fitting results of 14 Chandra GCS sources (F = 2.5 ) 



Obs.ID 


Source 


Counts 


iVif(10^^cm-2) 


xVdof 


2271 


J174722.9-280904 


108 


16.81 ±1.101 




16.926/21 


2273 


J174639.1-285351 


179 


(1.464 ± 0.378) X 10 


-4 


21.616/14 


2273 


J174622. 7-285218 


119 


15.02 ± 2.732 




18.401/11 


2274 


J174705.4-280859 


141 


0.2518 ±0.2823 




8.073/12 


2275 


J174319.4-291359 


133 


1.234 ±0.2266 




10.669/11 


2276 


J174550.4-284921 


180 


5.534 ± 0.641 




29.323/19 


2276 


J174550.4-284911 


86 


4.989 ± 1.992 




6.515/5 


2277 


J174804.9-282917 


94 


14.92 ± 3.407 




4.583/6 


2277 


J174729.0-283516 


204 


0.8724 ± 0.2423 




23.174/14 


2278 


J174417.3-293944 


695 


0.7984 ± 0.0907 




26.578/24 


2282 


J174602.2-291039 


174 


1.759 ± 0.252 




18.668/16 


2285 


J174626.1-282530 


158 


4.640 ± 0.616 




13.99/13 


2294 


J174502.8-282505 


107 


0.621 ± 0.250 




12.533/11 


2295 


J174451.7-285309 


104 


3.703 ± 0.655 




6.66/8 
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Fig. 3 Relation between Nu and HR' 
& HR2' . The solid line with squares is 
the relation of Nh and HR' , while the 
dotted line with triangles is the one of 
Nh and HR2' . 




Fig. 4 The Nh distribution of Chandra 
GCS sources. The black squares are the 
13 sources used for determining their 
joint photon index (§2, one of the 14 
sources cannot get the Nh value us- 
ing our method); their vertical locations 
have no meaning. Since their Nh values 
span across almost the whole range, it is 
believed that there is no selection effect 
involved. 



Assuming that F = 2.5 for all sources, Nh value for each source can be estimated 
according to the hardness ratios of the source. Following the definitions of Wang (private 
communication): HR = jq^, HR2 = where a, b and c are counts ^ in three 

bands: A (1-3 keV), B (3-5 keV) and C (5-8 keV), respectively. The values of HR 
and HR2 for each source have been given in the source list. Using the PIMMS software 
http: / /cxc. harvard.edu /toolkit /pimms.jspl) , the relation between HI column density Nh 
and hardness ratios can be obtained as shown in Fig. 3. Here we should emphasize that 
the definitions of hardness ratios in Fig. 3, labelled as HR' and HR2' there, are different 
from the ones above, HR and HR2. This is because when using PIMMS, we can only 
get the values of count rates, i.e. the parameters of Poisson processes in three bands 
Aa, As, Ac: Their definitions are HR' = k^-^^, HR2' = 4^^^. 

We therefore use the method in Jin et al. (2006) to estimate the values and error 
intervals of HR' and HR2' . For a Poisson process, the A parameter obeys the Gamma 
distribution under certain counts as follows, 

p{\A=x\nA=a) = ^^-^j — . (1) 
a! 

^ An off-angle efficiency correction has been applied. However the maximum correction factor 
for the 8.638 keV photons at 8' off-axis angle is only 1.67. Therefore the corrected counts still 
follow Poisson distribution closely. 
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Angular distance to the Galactic center (deg.) Angular distance to the Galactic center (deg.) 



Fig. 5 Angular distribution of Chandra GCS sources. Left panel: source count 
in each interval of angular distance [0°,0.05°), [0.05°, 0.10°) ... [1.20°, 1.25°); 
Right panel: source number density of each interval of angular distance 
[0°,0.05°), [0.05°, 0.10°) ... [0.95°, 1.00°). 



Based on this result, the probability density function of Xa/^b is, 

Aa , ^"(g + ^ + l)! 

pi-=z\nA = a,ns = b) = ^^y^,^,^^, - (2) 



The expectation of Xa/Xb is 



E{^\nA = a,nB = b) = ^4~- (3) 

Ab 

The most probable value zq can be obtained as, 

= bT-2- 

Hence, for a Poisson process, Xa/Xb is not exactly equal to a/b. Here the most probable 
value is taken as the estimate of Xa/Xb- For the hardness ratio defined as HR' = Xb+\a ' 
the probability density function is, 

.Xb-Xa I (l-z)"(l + z)''(a + 6+l)! 



Again the most probable value zo is taken as the estimate of 



Xb-Xa^ b-a 

^0 = ^x^Tx^^^'^^^ = bT^- 

The method of error estimation is shown as the following equation (the error interval 
denoted as [C, D]), 



JC Xb + Xa J-1 Xb + Xa 



(7) 



° Xb — Xa s, ,Xb — Xa 



Pif^—Y^ = ^)dz = 90% X / Pi^^—^ = z)dz. (8) 

iza Ab -\- Aa Jza + Aa 

According to the relations in Fig. 3, the Nh values and errors for each source can 
be derived by linear interpolation. As can be seen from Fig. 3, HR' and HR2' estimates 
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Fig. 6 The correlation of unabsorbed fliix of X-ray sources with HI column 
density. The bold line shows the detection threshold of this survey. The five 
crosses stand for the scales of the error bars of Nh and fiux values. The extended 
dotted lines show the Nh error intervals if the source photon indices T vary in 
2 ~ 3 range. 

are more accurate for Nh < 2 x 10^^ cm~^ and Nh > 2 x 10^^ cm~^, respectively; 
thus the Nh value for each source are assigned accordingly. Because the curves in Fig. 
3 cannot cover the whole range of HR' and HR2' as [—1,1], the Nh values of nearly 
80% 600) of the GCS sources are acquired in this procedure, which are used in the 
following investigation. The Nh distributions derived from HR' and HR2' are shown in 
Fig. 4. 

3 SPATIAL DISTRIBUTION OF GCS SOURCES 

3.1 Angular distribution 

Fig. 5 shows the angular distribution of these sources. For angular distance in the 
[0.25°, 0.30°) interval, there are the largest number of sources. However on the source 
number density curve, the Galactic center is the densest region. A sharp decrease exists 
in the immediate outer range. For the region with angular distance > 0.1°, the decrease 
of source number density is slower. The total decrease is slower than the exponential 
trend. 

3.2 Radial distribution — Where are the sources? 

Our goal is to find out the three dimensional distribution of these sources from the 
projected two dimensional image. It is natural to consider that the sources arc likely 
concentrated in the Galactic center region since that region is known to be crowded by 
various kinds of objects. Another possibility is that the major part of the sources are 
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distributed throughout the Galactic disk. Here Monte Carlo simulation is used to study 
that which model is consistent with the statistical properties of these sources. 

We start with the Fig. 6. The horizontal axis of this figure is the HI column density 
of these sources, while the vertical axis is the unabsorbed flux in 0.2-10.0 keV (all the 
fluxes hereafter are in the same energy band), calculated from the T = 2.5 power-law 
model, the corresponding HI column density and the count rate of each source. The blank 
region in the bottom right corner reflects the sensitivity of this survey. The criterion of 
the source detection is based on Poisson probability (Wang 2004). The probability of a 
count deviation above the background Cbg can be formulated as 



where Cu and Cbg are the total counts and background counts for each source, respec- 
tively. If the probability P is less than a preset threshold Pth, the source detection could 
be declared. Here in this Chandra GCS, the Pth is 10~^. In our case, this threshold could 
be approximated as the S/N > 2, where S/N = {Cu — Cbg)/ VCu- This point can be 
seen in the following simulations. The bold curve in Fig. 6 shows the detection threshold 
in this approximation. 

We design some simulated observations. We assume a set of sources with certain 
luminosity function, and have them spread on the disk or concentrated at the Galactic 
center. The parameters in the simulation are: 1) cumulative luminosity function: power- 
law with r = 0.5, which is similar to the discrete X-ray sources in nearby spiral galaxies 
(e.g. Colbert et al. 2004; Tennant et al. 2001); 2) HI column density distribution: For the 
case that these sources are concentrated in the Galactic center, we simply take a uniform 
distribution in order to minimize the number of simulation parameters; For the case that 
sources are spread throughout the Galactic disk, the HI column density is assumed to be 
proportional to the radial distance, taking the HI column density of the Galactic center 
as 1 X 10^^ cm~^ corresponding to 8.0 kpc distance (Baganoff et al. 2003); the farthest 
boundary of the Milky Way to us is within 30.0 kpc; 3) background count distribution: 
the same as the Chandra GCS observation, i.e., a Poisson distribution. 

Fig. 7 and Fig. 8 summarize the simulation results, test is used to find out which 
simulation is closer to the observation. For two discrete data sets, data in each set are 
divided into several intervals. The number of intervals is v. Rj and Si are the numbers 
of data in the zth interval for each set, respectively. Then the value is defined as the 
following equations (Press et al. 1992), 




(9) 




Ri + Si 



(10) 





(11) 
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Fig. 7 HI column density vs unabsorbed fliix, where bold curves in all three 
panels show the detection threshold of the real survey. Left panel: simulation by 
assuming sources are uniformly spread all over the Galactic disk. The value of 
HI column density is proportional to its radial distance. Middle panel: in real 
Chandra survey. Right panel: simulation by assuming most sources concentrated 
in the Galactic center. The value of HI column density is uniformly distributed. 
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Fig. 8 Left panel: Comparison of Nh value distributions in real observation 
and simulations. The simulation in source coneentrating scenario is closer to 
the observation. The main difference shows in low Nh value interval. Right 
Panel: Comparison of unabsorbed flux value distributions in real observation and 
simulations. Again the simulation in source concentrating scenario is more likely. 
The first point of each line stands for sources with unabsorbed flux < 10"^"^ ergs 
cm~^ s~^, while the last point stands for sources with unabsorbed flux > 10~^^ 



ergs cm ^ s ^. 



The degree of freedom for this distribution is v. 

First, we examine the case that X-ray sources are spread out in the disk. Fig. 8 shows 
that the simulated HI column density distribution and the unabsorbed flux distribution 
in this case are very different from those of the observation: ^ jdof = 199/40 for the Nh 
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distribution, ~)(^ /dof = 126/22 for the unabsorbed flux distribution. The main disagree- 
ment is that there are too many sources with HI column density at 2x 10^'^— 4x 10^^ cm~^ 
in this simulation than in observation. The reason is that the source number density does 
not vary with the Galactic latitude, i.e. the number density is uniform in the field of view. 
However the density of the Galactic material is known to decrease for high latitude re- 
gion. An exponentially decreasing factor exp(— 2/2:0) is then used. For smaller zq value, 
the /dof value is found to decrease: even for zo ~ 1 pc, the simulated distribution of 
HI column density cannot satisfactorily match the observation (x^/^o/ = 111/40). Such 
small value of is inconsistent with our knowledge of Galactic mass distribution: for var- 
ious kinds of matter, zq ^ 70 — 400 pc (Ferricrc 2001). The mass model of the Milky Way 
(Bahcall & Soneira 1980; Dehnen & Binncy 1998; Grimm, Gilfanov & Sunyaev 2002)is 
also tested and simulation cannot match the observation either [y^ /dof = 294/40 for 
the Nh distribution). Furthermore, there are many differences between the left panel 
and middle panel of Fig. 7. A void appears in the top right part of the simulation figure, 
which is not present in the observation. The reason for this void is that for the same 
luminosity range, the flux of a farther source is smaller. There is no reason to believe 
that farther sources are more luminous. Therefore this difference cannot be accounted 
for easily. In fact, we find that this kind of void exists whether or not we add the expo- 
nentially decreasing factor. We thus conclude that the sources in Chandra GCS are not 
likely distributed throughout the Galactic disk. 

Then we examine the other case: Point-like sources are concentrated in the Galactic 
center. The HI column density distribution is assumed to be a uniform distribution. 
Comparisons with the observed HI column density and flux distributions are: /dof = 
104/40 and ^ jdo^ = 83/22, indicating significant improvement over the case that these 
point sources are distributed throughout the galactic disk. However even this case cannot 
match the real observation completely. The middle and right panels of Fig. 7 show the 
differences. One major difference is on the top left part of the two panels. In observation, 
there is a void on the top left part, i.e. less high flux sources with small HI column density 
are detected than in simulation. This void cannot be explained by selection effect, since 
selection effect only makes soiirces under the threshold curve; undetectable. Therefore, 
from this point we infer that some kind of positive correlation exists between the HI 
column density and the unabsorbed ffux, i.e. brighter sources have higher HI column 
density. Since these sources are probably concentrated in the Galactic center, they should 
have approximately the same distances. We therefore conclude that a positive correlation 
exists between the X-ray luminosities and HI column densities of the Galactic center 
sources. The lack of high flux sources with low HI column density values might also be 
due to the pileup effect for low Nh sources or the cumulative luminosity function steeper 
than 0.5. However, since there are very few sources (less than 1%) which have significant 
pileup effect, this kind of void is not likely due to the pileup effect. 
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4 A POTENTIAL INTERPRETATION FOR THE NATURE OF 
THESE SOURCES 

As shown in the previous section, the luminosity of the X-ray sources is positively corre- 
lated to the HI column density. One possibility is that the X-ray luminosity comes from 
the interaction between an isolated neutron star and interstellar medium (ISM), mainly 
dense molecular clouds in the Galactic center region. 

Since stars in the Galactic center are very old, there should be many dead isolated 
neutron stars. They may move into the molecular clouds with speed ~ 10^ km s~^ 
(Arzoumanian et al. 2002), and then accrete the ISM via standard Bondi accretion. The 
statistical properties of molecular clouds in Galactic center region (Miyazaki & Tsuboi, 
2000; Oka et al, 2001) are as the following: 1) mass M: 1.3 x 10^ — 1.3 x 10* M©; 2) 
size R: 1.3 — 50 pc; 3) mass density p: 10~^^ — 10~^* g cm~^; 4) number density n: 
10^ — 10^ cm~^. The sound speed of the molecular clouds (~ 1 — 10 km s~^) is much 
smaller than the velocity of the moving neutron star. The accretion rate of standard 
Bondi accretion should be calculated as 

47r(GM„)^mpn 

(F2 + ^2^3/2 

where V7 is the velocity of neutron star taking unit as lO'' cm s~^. Here we take V7 = 1, 
neutron star mass M„ = 1.4 Mq, neutron star radius i?„ = 10 km, then the accretion 
luminosity should be 

iacc~^^~10^^nergss-^ (13) 

Substitute the number density of molecular clouds into the above equation, the luminosity 
should be lO"*^ — 10'^''' ergs s^^. For Chandra GCS sources studied above, the unabsorbed 
fluxes stay in the range of 10^^^ — 10^^^ ergs cm^^ s^^. If most of these sources are at 
the Galactic center with a distance of 8.0 kpe, the luminosities should be lO"*^ — 10'^''' 
ergs s^^, consistent with our result. The valuers of c;olumn density from the center to the 
border of the dense molecular clouds are 10^^ — 10^^ cm^^. Thus for the neutron stars in 
the molecular clouds, the HI column density should be in this range. Since the molecular 
clouds have different volumes, the HI column density is not strictly proportional to the 
number density of a molecular cloud. However, a positive correlation does exist. For 
sources with small HI column density, they are not likely to have high accreting X- 
ray luminosity. Therefore in this scenario, the luminosity of X-ray sources is positively 
correlated with the HI column density as in the observation. The top left void in the 
middle panel of Fig. 7 can be explained in this scenario. 

From the discussion above, we also argue that the absorption of X-ray is mainly caused 
by the dense molecular clouds near the Galactic center, rather than by the Galactic disk 
ISM along the line of sight to us. In fact, using the ISM model of the Milky Way (Ferriere 
1998, Eqs.(6)(7)) the HI column density from the Galactic center to us is calculated as 



12 J. F. Wu, S. N. Zhang, F. J. Lu & Y. K. Jin 

only ~ 10^^ cm~^, only a small fraction of the HI column density of Chandra GCS sources. 
Our model can also explain why the HI column density of sources in the Galactic center 
varies in a wide range. The HI column density of foreground sources should be < 10^^ 
cm~^. Wc can see in Fig. 6 that the foreground sources are only a small part of the 
whole catalog. The percentage of background extragalactic sources is also very small 
(Bandyopadhyay et al. 2005). Thus the contribution from foreground and background 
sources can be ignored. 

Treves et al. (2000) reviewed the theory and the observations of old isolated neutron 
stars in the Milky Way. They estimated that there should be 10^ — 10^ neutron stars 
in our Galaxy, and most of them are dead. The Galactic center region is thus a good 
place to observe old isolated neutron stars. Neutron stars in this region are so old that 
the rotation frequency and magnetic field have decayed and the neutron stars could be 
accretors, instead of ejectors and propellers. Even in ejector or propeller phase, accretion 
could also occur (Zhang, Yu & Zhang 1998): A small partion of materials accretes onto 
the polar region of the neutron star through a quasi-spherical ADAF. This scenario may 
occur in neutron star soft X-ray transients. For our model, the accretion onto isolated 
neutron stars can also produce variability due to the instability of magnetic field or the 
fine structure of ISM (Treves, private communication). 

For those neutron stars still alive (i.e. radio pulsars), there can also exist the positive 
correlation between their X-ray luminosity and the ISM density. As these neutron stars 
move into ISM with a speed much higher than the ISM sound speed, there should be a 
bow shock. Romani, Cordes & Yadigaroglu (1997) discussed the Guitar nebula formed 
in this scenario. The shock luminosity from the shocked ISM along the Guitar body is 
proportional to the hydrogen number density. Recent research on bow shocks also showed 
that for some region of the shock, the X-ray luminosity is positively correlated with the 
ISM density (Gaensler et al. 2004). 

Further study on the nature of point-like X-ray sources needs miiltiwavelength obser- 
vation. Bandyopadhyay et al. (2005) took a VLT infrared survey on 77 X-ray sources of 
the Chandra GCS field. They detected candidate counterparts for 75% of the sources in 
their sample, and suggested that the X-ray sources might be wind-fed accreting neutron 
star binaries. However, their sample is only a small part of the sources. Laycock ct al. 
(2005) of ChaMPlane program excluded High Mass X-ray Binaries (HXMBs) as the main 
contributor of X-ray sources in the Galactic center. Other kinds of sources account for 
more than 90% of the X-ray sources. Furthermore, Belczynski & Taam (2004) argued 
that accreting neutron star systems are not likely to be the main contributor to the faint 
X-ray sources in the Galactic center, because neither wind-fed nor Roche lobe accreting 
neutron star systems can explain the amount or spectral properties of the sources in 
catalogs of surveys in Wang, Gotthelf & Lang (2002). Our model of the X-ray sources. 
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old isolated neutron stars with Bondi accretion, can explain the faint soft X-ray source 
population in the Galactic center. 



5 SUMMARY 

We have performed statistical analysis to 600 of 800 point-like sources detected in 
the Galactic center survey of Chandra X-ray Observatory due to its excellent spatial 
resolution and broad energy range. Fourteen bright sources detected are used to fit jointly 
an absorbed power-law model, from which the power-law photon index is determined to 
be ~2.5. Assuming that all other sources have the same power-law form, we use the 
relation between hardness ratio and HI column density Nh to estimate the Nh values 
for all sources. We study the spatial distribution of these sources with Monte Carlo 
simulations and infer that these sources are mainly concentrated in the Galactic center 
region. The luminosities of X-ray sources are found to be positively correlated with the 
HI column density. We propose a model of old isolated neutron star with Bondi accretion 
as the X-ray sources to explain the correlation. This model is consistent with the present 
knowledge of dense molecular clouds in the Galactic center. We also argue that the 
molecular clouds should be the main contributor to X-ray absorption, rather than the 
ISM in the Galactic disk. Further extensive multiwavelength surveys are needed to study 
the nature of X-ray sources in the Galactic center further. 
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